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Spin and valley-orbit splittings are calculated in SiGe/Si/SiGe quantum wells (QWs) by using the 
tight-binding approach. In accordance with the symmetry considerations an existence of spin split- 
ting of electronic states in perfect QWs with an odd number of Si atomic planes is microscopically 
demonstrated. The spin splitting oscillates with QW width and these oscillations related to the 
inter-valley reflection of an electron wave from the interfaces. It is shown that the splittings under 
study can efficiently be described by an extended envelope-function approach taking into account 
the spin- and valley-dependent interface mixing. The obtained results provide a theoretical base to 
the experimentally observed electron spin relaxation times in SiGe/Si/SiGe QWs. 

PACS numbers: 73.21. Fg, 78.67.Do 

I. INTRODUCTION 

At present various semiconductor materials are being involved in the spintronics activities. SiGe/Si quantum well 
(QW) structures are among them. Silicon-based systems can be particularly promising due to a comparatively weak 
spin-orbit interaction and long electron spin-relaxation times. Although bulk Si and Ge have an inversion center, QW 
structures grown from these materials can lack such a center and allow the spin splitting of the electronic subbands, 
even in the absence of structure inversion asymmetryi An ideal SiGe/Si/SiGe QW structure with an odd number of 
Si atomic planes is characterized by the T>2d point-group symmetry and, therefore, allows spin-dependent linear-in-fc 
terms in the electron effective Hamiltonian 



H (1) (fc||) = a(cr x k x - (Tyky) 



(1) 



where <r x , a y are the spin Pauli matrices, fc|| is the two-dimensional wave vector with the in-plane components k x , k y , 
and x || [100], y || [010]. 

In the present work we use both the microscopic tight-binding model and the envelope- function approach to calculate 
the spin splitting of the conduction subbands in diamond-lattice QWs. The obtained results are of particular interest in 
connection with the experimental studies of electron spin relaxation in Si/SiGe heterostructures^^ The consideration 
of a Si/SiGe structure with perfect interfaces and without built-in electric fields allows one to put the upper limit to 
the electron spin relaxation time. 
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FIG. 1: Schematic representation of the lower conduction bands Ai (solid curve) and A2 (dashed) in bulk Si along the Y— X 
direction in the first Brillouin zone. Horizontal bars in figure illustrate extension of the el quantum-confined state in the k 
space; fci, k2, k[ and k' 2 are wave vectors of four Bloch states mixed in a QW. 
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In a bulk homogeneous sample of Si, two of six equivalent minima of the conduction band Ai are located in two 
points, k and — fc , along the direction [001] of the first Brillouin zone as illustrated in Fig. ^ The point-group 
symmetry of a Si/SiGc(001) QW reduced and allows mixing between four bulk Bloch states attached to the k and 
—kg valleysA^ The valley-orbit mixing occurs under electron reflection from a heterointerface: an electron with the 
wave vector fci sa fco is reflected not only to the state k 2 attached to the same valley fco but also to the state k' 2 in the 
second valley —fco, see Fig.^ The reflected wave is a superposition of two waves with their phase difference dependent 
on the distance z from the interface as 2koZ. In the QW grown along the [001] direction, quantum-confined electron 
states are standing waves formed as a result of multiple reflection of the four waves fci, k2, k[, k' 2 , or ±fco ± (fci — fco), 
from the both heterointerfaces. 

The spin splitting in conduction subbands is directly related to spin dependence of the electron oblique-incidence 
reflection from an interface. Spin-dependent reflection of an electron wave from interface consists of intra- and 
intervalley contributions. The latter should oscillate with the QW width L in the same way as the spin-independent 
valley-orbit splitting. Thus, interface-induced spin splitting A sp i n contains two contributions one oscillating with L 
and another being smooth. Their relation can be obtained in microscopic evaluations. 

The paper is organized as follows. In Sec. II we extend the envelope function method to take into account intra- 
and inter-valley spin-dependent contributions to the effective interface potential. In Sec. Ill we develop the sp 3 s* 
tight-binding model in order to calculate the dependence of the coefficient a in Eq. on the QW width, discuss the 
results of calculations and compare them with the analytical equations derived in Sec. II. The paper is concluded by 
Sec. IV. 



II. EXTENDED ENVELOPE FUNCTION METHOD 



Let us consider a QW layer A sandwiched between barriers B and C on the right- and left-hand sides, respectively. 
We assume that the three bulk materials j = A, B, C have the diamond-like lattice, the structure is grown along 
the principal crystallographic axis z \\ [001], and the lowest conduction subband el is formed by electronic states 
in the two A valleys with the extremum points ±fcoj = (0, 0, ±fcoj). Note that, in the Sii-^Ge^ solid solution, the 
extremum-point position is a function of the content x and values of koj are layer dependent. Because of the lattice 
constant mismatch some of the structure layers are strained. The layers B and C are assumed to be thick enough for 
the tunnelling tails of the quantum-confined el states to decay within these layers so that they can be considered as 
semi-infinite. 

In the generalized envelope function approximation the electron wave function 'J'(r) inside the layer j is written as 

*(r) = e ik rP[ ( p 1 (z;j)^ hoj (r) + cp 2 (z;j)^ ko] (r)} . (2) 

Here 

^o J (^)=e i ^ z u feoj (r) and ^ haj (r) = fT^'u^, (r) (3) 

are the scalar Bloch functions at the two A extremum points, u±k 0j (r) are the Bloch periodic amplitudes, (fi(z;j) 
and (f2{z;j) are the smooth spinor envelope functions defined within the layer j, p is the in-plane component of the 
three-dimensional radius- vector r. 

The two- valley effective Hamiltonian TL is presented as a sum of the zero-order valley- and spin-independent term 



H = - 



1 d k^, -\- ky 



(4) 



dz mi(z) dz m t (z) 

and the interface-induced ^-functional perturbation 

W = V L S(z -z L ) + V R S(z - z R ) . (5) 

Here mi and mt are the longitudinal and transverse effective masses for electrons in the A valley, zi and zr are the 
coordinates of the left- and right-hand side interfaces, Vl and Vr are both valley- and spin-dependent operators. Here- 
after we assume that the latter contain no differentiation d/dz, this assumption excludes the need in symmctrization 
of Vl,r and the 5-function. 

The form of Vl, Vr can be specified by applying the symmetry considerations. A single (001) interface is charac- 
terized by the C2v point-group symmetry allowing two linear-fc|| spin-dependent invariants, namely, 

h(k) = a x k x — (Tyky and h'(k) = a x k v — a v k x . 
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It follows then that the matrices V m (m = L, R) acting on the bispinor vector (pi.1/2, ¥1,-1/2, ¥2,1/2, V2.-1/2) can be 
presented in the form of a 2x2 block matrix 



V m = 



S rn h{k) + S' m h'{k) 

A; n i + p; n h(k) + p£h>(k) 



A m I + P m h{k) + P^h'(k) 
S m h{k) + S' m ti '(fe) 



(6) 



with its components being linear combinations of the Pauli matrices and the 2x2 unit matrix /. Here k = k 
and S m ,S' m , A n ,P m ,P m arc coefficients characterizing the right-hand (m = R) and left-hand (m = 
the first two of them (S m , S' m ) are real while others are complex. The diagonal components V m -n 



L) 



I- 



interfaces, 
V m -22 give 



Kn-2i describe interface-induced inter-vallcy 



intra- valley contributions whereas the off-diagonal components V m: r 

mixing. It is more convenient to perform the further considerations for a particular case of coinciding barriers, C = 
B, and coinciding extremum points, feoB = &oa (or fcoB = ^oa)- Then we briefly discuss how these considerations are 
generalized with allowance for C ^ B and different positions of extremum points koj . 

The choice of the electron Hamiltonian in the form of Eqs. @, J5j corresponds to a particular set of boundary 
conditions. For the structure B/A/B with fcoe = fcoA = fco; this set reads 

ip(z L + 0) = ip(z L - 0) , ip(z R + 0) = tp(z R - 0) , 



1 (dip 
m/(B) \ dz 



1 (dip 
mj(A) \dz 



ZL+O 



2 
n 



(7) 



1 (dp) 
mj(B) Vdl 



1 /d^ 

mj(A) ydz 



2 ^fl <P{zr) 1 



where ^(zi,^ ±0), [dip / dz) ZL r ±q are the envelope function and its first derivative at z approaching the interface L, R 
from the right- (+0) and left-hand (—0) sides. 

The next step is to analyze the phases of the coefficients A m , P m , P' n in the off-diagonal components of V m and 
establish a relation between Vl and Vr. First of all, we take into account that the translation of the radius- vector r 
by a three-dimensional Bravais-lattice vector a, results in a multiplication of the Bloch functions 4 ! ±k oj (f) in Eq. (0) 
by the factors exp (±ikoja z ), respectively. Therefore, one can present the coefficients in the off-diagonal components 
of V m a*££ 



A„, — X r 



-2iknz n 



3 P m Pm^ 



-2iknz n 



-2ikoz rl 



(8) 



where the complex coefficients \ m ,Pm-,p' m are independent of the interface position. In the following we assume the 
origin z = to lie in the QW center. 

The structure B/A/B is invariant under the mirror rotation operation £4 with the transformation center at z = 0, 
if the number N of atomic planes in the layer A is odd, and under the space inversion operation i, if N is eveni This 
symmetry property allows one to establish the relations between the coefficients in Eq. © for the left- and right-hand 
side interfaces. Since each of the operations results in the reciprocal transformation ipk oj (r) <-> ip-k oj (r) one has 



" A L " 




'0 1" 




' A R " 




'0 1" 


X* L 




1 




. A R . 




1 



or, equivalcntly, Al = Ar. Taking into account that, under the mirror-rotation operation £4, the C2t,-group invariants 
h(k) and h!(k) transform, respectively, into h(k) and —h'(k) while, under the space inversion i, both h(k) and h'(k) 
change their sign, we also obtain the relations 



and 



Sr = Sl , S'r = -Si , pr = pi , p' R 



Sr = -Sl , S R = -Si , p R = -pi , p' K 



-Pl 



-P'i 



(odd N) 



(even N) 



(9) 



(10) 



Hereafter we use the notations A, S, S' ,p,p' instead of Xr, Sr, S R ,pR,p' R . By using Eqs. © and HlOfl we can reduce 
the components in the matrix © to 



V R , 



Vr, 22 = Sh(k) + S'h'(k) , V U11 = VL l2 a = Sh(k) - S'h'(k) , 



(11) 
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FIG. 2: Schematic representation of (a) hierarchy of the el subband splittings and (b) the Si/Sii-ajGea; structure under 
consideration. We remind that, in this structure, the conduction band offset is mostly determined by the strain; 



Vk,i2 = Vl 21 = c- ik « L [\I+p h{k)+p'ti{k)} 



Vl.u = VS 21 = e lkoL [\*I + p*h(k)~p'*h'(k)} 



if N is odd and to 



Vk,ll = V R ,22 = Sh(k) + S'ti{k) , V^xi = V L , 



-Sh(k) - S'h'(k) , 



(12) 



Vr,i2 = Vl 21 = c~ lk " L [XI+p h(k)+p'h'(k)] 



Vl,i2 = K,2i = e lk ° L [\*I-p*h(k) -p'*ti(k)] 

if N is even. Here L — zr — zl is the QW width, it is given by L = Na /4 with being the zinc-blende lattice 
constant. 

Equations i|ll|) and (|12f> present the results of the extended envelope-function method and yield relations between 
coefficients in the matrices Vl and Vr for macroscopically symmetric QWs. 

If the barriers are grown from different materials B and C then the coefficients in Eqs. Qll[l. (|12|l should be labeled 
by the interface index, C/A or B/A, e.g., S(C/A) and p'(B/A). The different positions of the extremum points koj 
are easily taken into account by replacing ipi(z;j) and (f2(z;j) (j = B, C) in Eq. J5J) and in the boundary conditions 
©by 

Vi(z;j) = e^- k ^ Vl (z;j) , <p 2 (z;j) = c-^-^^z; j) , 

where Zj is the coordinate of the interface between the layers A and j = C or B. This replacement allows to retain 
the form of the perturbation TL' defined by Eqs. JSJ, (JSJ and (jSJ. 



A. Valley-orbit splitting 

The numerical calculations presented in the following sections confirm the hierarchy 

E x - E(k ) > E el > A v _ > A spin = a±k (13) 

illustrated by Fig. Here E(ko) and Ex are the conduction-band energies at the extremum point fco and the X 
point in the bulk material A, E e \ is the quantum-confinement energy for the lowest conduction subband, A v _ and 
A sp in are the valley-orbit and spin splitting of the el-subband states. Therefore, we can line up the discussion in series 
starting from the quantum confinement, turning then to the valley-orbit splitting and finally to the spin splitting. 
As above we start from the analysis of the symmetric structure B/A/B shown schematically in Fig. ^jp and then 
generalize the results on asymmetric structures with different barriers B and C. 

For eigenstates of the zero-approximation Hamiltonian 7io the inter-valley mixing is absent and the envelope func- 
tions referred to the first and second (OOl)-valleys form identical sets. In particular, for the el subband states in the 
B/A/B structure, the envelope has the standard form 

y(z)-c( C ° SqZ > if l z l- L / 2 ' (14) 

x[z) ~ c \ cos(q£/2)exp[-a8(|z| - L/2)] , if \z\>L/2. [ ' 

Here q = [2m l (A)E e i/h 2 ] 1 ^ 2 , ac = [2mi{Q){V - E^/h 2 ] 1 / 2 and c is the normalization factor. The size-quantization 
energy E e i satisfies the transcendental equation tan(gi/2) = (?e/q)[mi(A)/mi(B)]. 
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Now wc switch on the inter-valley mixing taking into account zero-fc terms in Eqs. I|llf) and (|12fl proportional to A 
and A* . According to Eqs. © , and l|12fl the matrix element of the inter- valley coupling is given by 

M M;2/ = |x (i/2)| 2 (Vr,i2 + V h , 12 ) = 2 |x (L/2)\ 2 |A| cos (k L - tj> x )8 at , , (15) 

where |A| and (f>\ are the modulus and the phase of A, and s,s' = ±1/2 are the electron spin indices. Thus, the 
energies of the split el states at k x = k y = are 

E el>± = ±2 \ x (L/2)\ 2 ■ |Acos (k L - <f> x )\ , (16) 

and the envelopes are 

¥>i(z; el, ±) = ±W2{z; el, ±) = x{z)/V2 , (17) 

where x{ z ) is defined in Eq. ifH)l and ry = signjeos (k L — </>>)}■ Therefore, the parity of the lower state |el, — ) (with 
respect to the operation 1S4 if N is odd and i if N is even) follows the sign of 77 and reverses with the reversal of 
77. Equation 1)16(1 expresses an oscillating character of A v _ Q (i) in terms of the envelope function method. It will be 
shown in Sec. Ill that this agrees with the tight-binding numerical results. 

For an asymmetric C/A/B structure, the inter- valley matrix element and el-subband energies are generalized to 

M M;2>S , = |x (L/2)\ 2 (\ R e- ik ^ L + X L e ikoAL ) 5 SS , , (18) 

E eh± =±\ X (L/2)\ 2 [ |A R | 2 + |A L | 2 + 2|A R A L | cos {2k M L + Al - Ar ) }^ 2 , 
where 4>\ m is the phase of A TO (m = R, L). 



B. Spin-orbit splitting 

The next step is to take into account spin-dependent terms in V m . Since the symmetry forbids spin-splitting of the 
electron states in the B/A/B system with an even number of atomic planes in the A layer, we set N to be odd. Then 
the inter-valley mixing is described by the matrix elements 

M M;2 , S , = 2 |x (L/2)\ 2 [ I A| cos (k L - cf> x )S ss , + \p\ cos (k L - <j> p ) h ss , (k) - i\p'\ sin (k L - <j> p , ) h' ss , (fc) ] , (19) 

where </> p , 4>p' are the phases of p and p' . Assuming the valley-orbit splitting to exceed the spin-orbit splitting we are 
able to rewrite the Hamiltonian in the basis (|17|) and obtain the following 2x2 spin-dependent effective Hamiltonians 
in the subbands (el, ±) 

H'(k; el, ±) = |x (L/2)\ 2 [V h , n + Vr,h ± V Re(F L , 12 + V K , 12 )} , (20) 

and finally 

H(k;el,±) =E elt± + a±h(k) , (21) 

where the coefficients in the linear-fc term are given by 

a± = 2 |x (L/2)\ 2 [ S ± \p\ V cos (k L - <j> p ) ] . (22) 

While deriving Eq. (|22|l we took into account both the intra- and inter-valley contributions to V m and retained only 
the terms up to the first order in S and p. In agreement with the symmetry arguments, neither the S"-dependent nor 
p'-dependent contributions to V m give rise to linear-fc terms. Note that, for the sake of completeness, in addition to 
the linear-fc terms one can include in the right-hand side of Eq. (|22|) a spin-independent quadratic-fc term h 2 k 2 /2m^. 
Here m^ 1 = (el|mjT (z)|el) and the angle brackets mean averaging over the el state defined in Eq. (|14|) . 

In addition to a smoothly decreasing term in a± predicted in Ref. 0, Eq. (|22|l contains an oscillating term. The 
reason for the oscillations is mixing of valley states at the QW interfaces. Tight-binding calculations presented below 
show that \p\ > S, i.e., the oscillating part of a± is dominating. 

For an asymmetric structure C/A/B, the linear-fc contribution to the Hamiltonian H(k; el,±) takes the form 

U (1) (fc; el, ±) = a ± h{k) + (3±h'{k) 

with 

a± - |x (L/2)\ 2 [ S R + S h ± Rc{c- i «(p R c- ifcoi +p L c ik « L )} } , (23) 

f} ± = |x (L/2)\ 2 [ + S[± Reie^Hp^^+p'^)} ] , 
and £ = arg {ARe _lfc ° L + ALC lfc ° L }. For the symmetric structure, e 1 ^ = J] = ±1. 
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III. TIGHT-BINDING CALCULATIONS AND DISCUSSION 

In order to estimate values of spin and valley-orbit splittings we have performed calculations of the electron disper- 
sion in the el conduction subband by using one of the empirical tight-binding models. More precisely, we have fixed 
on the nearest-neighbor sp 3 s* tight-binding model optimized for the conduction band^ This model is a reasonable 
compromise between the numerical load and the accuracy of representation of the band structure. It is capable to 
reproduce the indirect gap although shifts the position of the conduction band minimum from the experimentally 
measured point fco = 0.85 x 27r/ao to the point ko = 0.62 x 2ir/ao. Note that a value of kg is hardly reproduced 
even in the more sophisticated methods, namely, the second-nearest neighbor sp 3 s* and nearest neighbor sp 3 d 5 s* 
tight-binding models ,i2iAi leading to fcoao/27r = 0.758 and 0.813, respectively. The applicability of the sp 3 s* model is 
confirmed by the fact, see below, that values of the valley-orbit splitting A v _ calculated in this work and by using 
the sp 3 d 5 s* model& are of the same order of magnitude. 

The empirical sp 3 s* tight-binding method was previously applied for calculation of the spin splitting in bulk GaAs 
and GaAs-based QWs>i£ The linear-in-k splitting in a QW was compared with the cubic spin splitting in bulk GaAs 
where the component k z was replaced by ir/dcaAs with dc&As being the width of the GaAs layer. The agreement was 
obtained after replacing dc a As by an effective value d^ aAs and adjusting the coefficient 7 in the cubic-in-k contribution 
to the electron effective Hamiltonian Tt^ (k). The need in the introduction of the effective parameters c^g^as 

and 7 ett 

can be related to an additional contribution to Ti^(k) coming from the reduced symmetry of interfaces, or in other 
words from the anisotropic orientation of interface bondsii In contrast to the zinc-blcnde-lattice heterostructures, 
in diamond-lattice QWs the Hamiltonian Ti^(k) has no bulk inversion asymmetry term proportional to 7 and is 
contributed only by the interface inversion asymmetry term described by the coefficient a in Eq. (p~|) - 1 

In the tight-binding method the electron Hamiltonian is presented by a set of matrix elements taken between atomic 
orbitals. If a heterostructure is grown from diamond-like semiconductors along the [001] principal axis one can write 
the tight-binding free-electron wave function 

ip = X] c «^ $ ™^fc( r ) ( 24 ) 

in terms of planar orbitals 

$«,,,fe(r) = J2 c ik rm ®Ar - r m ) . (25) 

tn 

Here n — 0, ±1,±2... is the number of atomic planes perpendicular to the growth direction z || [001], § v is the 
orthogonalizcd atomic orbital with v being the orbital index, the index m enumerates atoms in the n-th atomic 
plane, r m is the position of the mth atom in this plane, in particular, z m = nao/4, k is the two-dimensional in- 
plane electron wave vector. The index v runs through 2N values where N is the number of orbitals taken into 
consideration and the factor is due to electron spin. For convenience we use below the Cartesian coordinate system 
x' || [110], y' || [110], z || [001]. In the nearest-neighbor approximation we obtain the following set of equations 

Ul,(2l)C 2 i-i + E (2l)C 2l + U x ,(2l)C 2l+1 = EC 2l , (26) 
Ul(2l - l)C 2 i- 2 + E {21 - l)C 2i _! + U y \2l - l)C 2l = EC 2l -i 

for the vectors C n containing 2N components c n ^. Here Eq{ji) is fe-independent diagonal matrices, U x i and U y > are 
k x i and k y i dependent matrices. The diamond lattice has two atoms per unit cell and can be represented as two 
face-centered cubic sublattices shifted with respect to each other by ^/Sao/i along the [111] direction. The atomic 
planes with even n = 21 and odd n = 21 + 1 (1 = 0, ±1...) belong to the different sublattices and differ in the direction 
of chemical bonds. As compared with the pair of planes 21 and 21 + 1, the orientation of chemical bonds between 
atoms in the planes 21 — 1 and 21 are rotated around the axis z by 90°. For brevity we omit here the detailed form 
of matrices U x ^ y i(n); for k x i = k y i = these matrices can be readily obtained from those for the zinc-blende-based 
heterostructures given in Ref. [13 ]. The matrices Eq, U x '. y ' are formed by the tight-binding parameters, which arc 
usually extracted from fitting bulk-material band structure to experimental one. The tight-binding parameters for 
Si and Ge are listed in Table [I] The diagonal energies are referred to the valence band top of each material. The 
parameters for Si were taken from Ref. those for Ge are not so critical for the purpose of this work, we collected 
them from Ref. ^5| and added a value of 0.30 eV for the spin-orbit splitting of the p orbitals ^lj. For SiGe alloys, 
we have used the virtual crystal approximation and the linear interpolation of the tight-binding parameters. The 
strain was taken into account only by shifting the diagonal energies E l/ in Si or Ge by the same value, the strain- 
induced splitting of the p-orbital states was ignored. The shift of diagonal energies for the barrier material is equal 
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TABLE I: Tight-binding parameters used in the calculations in eV. 





E s 


E p 


E s , 


Vss 




V xy 


v sp 


v s , P 


A 


Si 


-3.65866 


1.67889 


3.87576 


-7.97142 


1.69558 


23.32410 


8.87467 


5.41174 


0.045 


Ge 


-5.88 


1.61 


6.39 


-6.78 


1.61 


4.90 


5.4649 


5.2191 


0.30 
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FIG. 3: Valley-orbit splitting A v - in Sii-^Ge^/Si/Sii-zGej; (x — 0.25) QW versus the number of Si mono atomic layers. 
Analytical results shown by crosses calculated using Eq. JTH with |A| = 385 meV- A, <jf> A = 0.3tt. 

to AE C — A£" 9 , where AE g is the difference in the band gaps of the well and barrier bulk materials and AE C is the 
conduction-band offset. For a Sii-^Ge^/Si/Sii-^Gca; QW st ructu re with the strained Si layer and the Ge content 
x = 0.25 we used a value of AE C = 0.15 eV relying on Refs. 0I3HII3- 

Squares in Fig. [21 show results of tight-binding calculations of the valley-orbit splitting A v _ in symmetrical 
Sio.75Geo.25/Si/Sio.75Geo. 25 QWs as a function of the number N of Si atomic planes sandwiched between the thick 
barriers Sio.75Geo.25- The valley-orbit splitting exhibits pronounced oscillations with the increasing QW width, in 
agreement with Ref. The oscillation periods in Fig. of the present work and in Fig. 3 of Ref. @ vary consid- 

erably due to the difference in values of fco obtained in the sp 3 s* model used here and the sp 3 d 5 s* model. However, 
the splittings A v _ are of the same order of magnitude, e.g., at N ~ 60 the oscillation amplitudes differ only by a 
factor of ~2 which can be explained by the obvious sensitivity of A v _ to the model used. 

Crosses in Fig. [21 represent the calculation of A v _ in the envelope- function approximation, Eq. I|16|> . with fco = 
0.62 x 27r/ao- While calculating the electron envelope function at the interface, x(E/2), we used values of V = 150 meV 
for the conduction-band offset and of 0.907mo (mo is the free electron mass) for the longitudinal effective mass m;(A) 
as obtained in the sp 3 s* tight-binding model optimized for the conduction bandj^ and, for simplicity, took m/(B) 
equal to m;(A). The modulus |A| and the phase <p\ were considered in Eq. i|16|) as adjustable parameters. Their best 
fit values turned out to be |A| = 385 meV- A, 4>\ = 0.3n. It is seen from Fig. [3] that the simple analytical theory 
developed in Sec. II is in complete agreement with the results of more sophisticated tight-binding calculations. 

In Figs. the spin-orbit splitting for the two valley-orbit subbands .Eel,- an d E e i i+ are presented. This is the 
first calculation of the spin-splitting, no previous theoretical estimations are available in order to compare with. We 
define the splitting A sp j n in Fig. as the energy difference between the states with the spin parallel and antiparallcl 
to the x axis. Then if the antiparallel state lies higher the sign of A sp j n is negative as in case of the upper valley-orbit 
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FIG. 4: Spin-splitting of the valley-orbit split subbands in a Sio.7sGeo.2s/(Si)jv/Sio.75Geo.25 QW with N = 25 as a function of 
the in-plane wave vector for k || [100]. Curves 1 and 2 correspond to the subbands -B e i,+ and -E e i,-, respectively. The definition 
of the sign of A sp i n is given in the text. 



split subband E e \,+ , see curve 2 in Fig.0J The calculation shows that, up to k < 10 6 cm 1 , the linear dependence 

A spin (fc) = a±k (27) 

holds, in agreement with the Hamiltonian (|22|l . It is the variation of a± with odd N which is shown in Fig. [5] As one 
can see from to Figs. I2HE1 the valley-orbit and spin splittings are conveniently presented in the meV and [ieV scales 
confirming our assumption (|13fl . 

Figure [S] shows that the spin splitting A sp i n is an oscillating function of the QW width. This demonstrates that 
the inter- valley spin-dependent mixing at the interfaces prevails over the intra- valley contribution to a±. Squares and 
diamonds in Fig. show results of tight-binding calculation. The spin splitting is plotted only for odd number of 
Si monoatomic planes because, for even TV, A sp ; n in the symmetric structures vanishes. Conventional and x-shapcd 
crosses are obtained as the best fit using Eq. ilL'L'l) and choosing the same values for fco and <fi\ as in Fig. [3] and the 
additional adjustable parameters \p\ = 0.53 • 10 -5 cV-cm 2 , <f) p — 0.55, S — 0.15|p|. 

Now we compare the value of a_ estimated in this work with that extracted by Wilamowski ct alp 2 " from spin- 
resonance measurements in a Si/Sii_a;Ge x QW structure with x = 0.25. Note that the value a_ = 0.55 • 10~ 12 eV-A 
presented in this reference for a 120A-thick QW should be decreased by a factor of 1.6, i.e., in fact a_ = 0.34 • 10~ 12 
eV-A, see Ref. Il8|. Our estimation of a_ gives a value smaller by a factor ~ 6. This means that in the sample 
studied in Ref. |2| the Rashba (or structure-inversion asymmetrjiiS) contribution to the spin splitting dominates over 
the intrinsic contribution considered here. Nevertheless, the experimental value of the spin splitting is not so far from 
the limit for a perfect QW structure. 



IV. CONCLUSION 



The sp 3 s* tight-binding model has been developed in order to calculate the electron dispersion in heterostructures 
grown from multivalley semiconductors with the diamond lattice, particularly, in the Si/SiGe structures. The model 
allows one to estimate the orbit-valley and spin-orbit splittings of the electron quantum-confined states in the ground 
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QW width (A) 



27 41 54 68 81 95 109 




Number of atomic layers 



FIG. 5: Spin-splitting constant a in Eqs. 12 II . I|27|l versus the QW width determined by the number of Si monoatomic layers 
(odd N are taken in consideration only). The spin splitting of the lower subband E e \- is shown by diamonds (tight-binding 
calculation) and x-shaped crosses (envelope function approximation), those for the upper subband £/ e i+ are shown by squares 
and conventional crosses. 



subband. In the employed tight-binding model, the spin-orbit splitting is mostly determined by the spin-dependent 
orbit-valley mixing at the interfaces. For this reason the coefficients a± describing the lincar-in-fc splitting are strongly 
oscillating functions of the odd number, N, of the Si monoatomic layers. 

In addition to the numerical calculations, an envelope-function approximation has been extended to take account 
of spin-dependent reflection of an electronic wave at the interface and interface-induced intervalley mixing. The 
dependencies of the valley-orbit and spin-orbit splittings upon the number of Si atomic planes calculated in the tight- 
binding microscopic model are successfully reproduced by using simple analytical equations derived in the envelope- 
function theory and fitting the parameters that enter into these equations. It follows then that the envelope-function 
approach can be applied as well for the description of electron-subband splittings in a realistic Si/SiGe structure. 
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